Fermi liquid behaviour in an underdoped high T c superconductor 
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We use magnetic quantum oscillations in the underdoped high T c superconductor YBa2Cus06+cc 
(x ~ 0.56) measured over a broad range of temperatures 100 mK < T < 18 K to extract the 
form of the distribution function describing the low-lying quasiparticle excitations in high magnetic 
fields. Despite the proximity of YBa2Cus06+cc (x ~ 0.56) to a Mott insulating state, various 
broken symmetry ground states and/or states with different quasiparticle statistics, we find that 
our experimental results can be understood in terms of quasiparticle excitations obeying Fermi-Dirac 
statistics as in the Landau Fermi liquid theory. 

PACS numbers: 71.45.Lr, 71.20.Ps, 71.18.+y 



At the heart of our understanding of itinerant electron 
systems is Landau's Fermi liquid theory in which gapless 
excitations on a Fermi surface in momentum space are 
governed by the Fermi-Dirac distribution 0, 0, 0, S] • O ne 
of the foremost examples of a breakdown of the Fermi 
liquid paradigm is in fractional quantum Hall systems [5[ . 
The validity of Landau Fermi liquid theory could also be 
in question in other interacting systems [6[ such as those 
at the border of electron localisation. 

The functional form of the Bose-Einstein distribution 
has been demonstrated in detail by means of a number 
of experiments including recent optical measurements in 
cold atomic gases [7]. Interestingly, however, it is not 
clear as to what extent of detail the precise form of 
the Fermi-Dirac distribution function - which is the ba- 
sis of the theory of electron systems - has been mea- 
sured [4]. In this paper, we test the validity of Fermi- 
Dirac statistics governing the elementary excitations in 
one of the pre-eminent, most strongly correlated exam- 
ples of a condensed matter system: the underdoped high 
T c cuprates [8[. Quantum oscillations provide a unique 
means of accessing the statistical distribution governing 
quasiparticles of an interacting system. While other ther- 
modynamic quantities such as heat capacity are useful 
in accessing the second moment (i.e. the variance) of 
the statistical distribution governing the system of inter- 
acting particles, we show here the quantum oscillation 
measurements access not only the second order moment, 
but also the higher order moments, and hence provide 
a more complete test of the applicability of the Landau 
Fermi liquid description. 

We measure magnetic quantum oscillations in 
YBa2Cu306+x (x ~ 0.56) [1, [lO] over an extensive range 
of temperatures (100 mK < T < 18 K) to obtain the tem- 
perature dependence of the oscillation amplitude a(T). 
Quantum oscillations originate from changes in the occu- 
pancy of partially filled Landau levels at the Fermi energy 



£f as the magnetic induction B is varied, making their 
amplitude particularly sensitive to changes in the thermal 
distribution function f(z) (where z = [e — £f]Ab^) H. 
On increasing the temperature T, the broadened proba- 
bility distribution \f'(z)\ = —di/dz smears the phase of 
the quantum oscillations in energy e (see Fig. [I]), caus- 
ing a(T) to be reduced. The function a(T) measures the 
Fourier transform of the statistical probability distribu- 
tion |f'(z)|- We therefore use the inverse fourier trans- 
form of the measured a(T) to access the distribution of 
interacting quasiparticles in underdoped cuprates. 

An accurate determination of |f'(z)l requires a(T) to 
be measured over a wide range of T. Here, we achieve a 
dynamic range of ~ 50 dB (see data in Fig. [2]) by using 
the contactless conductivity technique (details described 
elsewhere pj], EH Ell ) in a motor generator-driven mag- 
net delivering slowly swept magnetic fields to ~ 55 T in 
4 He medium down to ~ 1 K and torque measurements 
in a continuous DC magnetic field reaching ~ 45 T in 
a dilution refrigerator. For measurements in the motor 
generator-driven magnet, the uncertainty in T is min- 
imized by ensuring that the sample is well-coupled to 
the liquid cryogen (i.e. immersed in liquid 4 He), thereby 
minimising heating due to irreversible field effects during 
the measurements for temperatures below 4.2 K. At these 
temperatures, the error in T is estimated from a compar- 
ison with DC field measurements, and with the resistive 
crossover field providing an accurate secondary indication 
of T. At temperatures above 4.2 K, the error in T is esti- 
mated by repeating field sweeps at different vapour pres- 
sures of 4 He gas. Figure [2] shows the measured quantum 
oscillations in resonant frequency shift [l2[ after poly- 
nomial background subtraction at various temperatures; 
performing a Fourier transformation of the oscillations in 
1/B using a Hann window yields a(T) of the most promi- 
nent (a) oscillation shown in Fig. OH Despite the close 
proximity of two other frequencies [14] , the temperature 
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FIG. 1: Considering the Fermi-Dirac distribution fFD = 
(1 + e z )~ 1 (where z = [e — £f]/&bT), the T-dependent step 
in occupation number (red/green lines) causes the oscillatory 
density of states g — goe l27re ^ huJc (assuming that go is approx- 
imately constant on the scale of the cyclotron energy %oo c = 
TieB/m*) shown by the blue line to be thermally smeared 
by the probability distribution |fpD(^)l — 1/2(1 + coshz) 
(red/green shaded regions). The consequent reduction in am- 
plitude is equivalent to a Fourier transform of |fpD(X)l? yield- 
ing oscillations oc e^ 27rF / s ^ periodic in 1/B (where F is the 
conventional quantum oscillation frequency |3|]) modulated 
by a T-dependent pre-factor a(T) = cionr// sinh ttt/ (where 
rj = 27vkBTm* /TieB and ao is a constant). The quantum os- 
cillatory magnetisation and resistivity of interest in this paper 
can be expressed in terms of the above thermally averaged 
density of states, and hence the same thermal amplitude fac- 
tor a(T). 
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FIG. 2: Magnetic quantum oscillations measured in 
YBa2Cu306+x with x ~ 0.56 (after background polynomial 
subtraction). This restricted interval in B = |B| furnishes a 
dynamic range of ~ 50 dB between T — 1 and 18 K. The 
actual T values are provided in Fig. [3] 



dependence of the a oscillation amplitude is unaffected 
since the other frequencies are much weaker than the a 
frequency, and the effective mass of all three frequencies 
has been measured to be nearly identical [141 ]. 

To extend the range in T to lower values, we include 
low temperature quantum oscillation data measured us- 
ing magnetic torque in a dilution fridge in DC fields. The 
temperature error in these dilution fridge measurements 
could arise either from a nonlinear thermometer magne- 
toresistance, or a drift in temperature stability during the 
magnetic field sweep. We calibrate the temperatures in 
a magnetic field by measuring the magnetoresistance be- 
tween 11 T and 45 T of the Ru02 thermometer used for 
these measurements, and collating them with the pub- 
lished values of magnetoresistance up to 8 T [15]. We 
find the magnetoresistance to behave very close to lin- 
earity for all T < 1 K, signalling that the chief error in 
T arises from thermal drift— estimated from measured 
changes between rising and falling field. Although the 
DC fields extend only up to 45 T (lower inset of Fig. [3]), 
the dependence of the amplitude on T/B is found to be 
independent of the interval in 1/B (upper inset of Fig.[3|), 
enabling the contactless conductivity and torque data to 
be combined (upon amplitude renormalization at 1 K) in 
a plot of the oscillation amplitude versus the dimension- 



less parameter r] = 2-Kk^Tm* /TieB in Fig. [4^. 

The quasiparticle probability distribution is obtained 
by performing a Fourier transform of the quantum oscil- 
lation amplitude: 



|f(*)i 



l (2tt) 



z a(T)dr] 



(1) 



Its comparison with the model Fermi-Dirac distribution 
|fp D (z)| involves two variable parameters: the effective 
mass ra* and the amplitude renormalisation factor ao- 
These are optimized by fitting to the even partial mo- 
ments 
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r(z)\dz 



(2) 



to those of the model probability distribution from the 
lowest order moment upward. While the second moment 
\±2 (i.e. variance) can be accessed by other thermody- 
namic experiments (e.g. heat capacity), a key advan- 
tage of quantum oscillation measurements is the abil- 
ity to compare the entire hierarchy of moments with 
probability distribution models. The chief experimen- 
tal limitation of quantum oscillation measurements is 
the breadth of the temperature range over which oscilla- 
tions are measured. To perform the Fourier transform, 
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FIG. 3: Quantum oscillation amplitude a versus T extracted 
by Fourier analysis of the data in Fig. [2] for 40 < B < 54 T, 
corresponding to an average B = 45.96 T. The upper inset 
shows the same data (filled circles) plotted versus T/B to- 
gether with similar points extracted from the field interval 
30 < B < 45 T in which the amplitude has been renormal- 
ized (open circles), matching the field range of the DC field 
torque measurements (lower inset). The close correspondence 
between field ranges implies that the dependence on T/B is 
independent of the interval in B within the resolution of the 
experiment, suggesting a magnetic field-independent m* . 



evenly spaced points are required in ±rj. These are ob- 
tained by symmetrizing the data in Fig. 0^ with respect 
to T = and making a linear interpolation (in the range 
— ?7iim < Vm < ?7iim ~ 1-8). The data is further padded 
with zeroes to yield a Fourier transform with points finely 
spaced in z. The combined effects of the finite spacing be- 
tween data points in Fig. [3] and experimental uncertainty 
restrict the reliable range in z to — z\i m < z < z\i m ~ 7.5— 
hence the comparison of 'partial' moments in Table 1 and 
Eqn.El 
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TABLE I: Evaluation of the lowest even moments. The first 
row contains the lowest even moments of the Fermi-Dirac 
probability (corresponding to zn m = oo). The second row 
contains the corresponding partial moments of the Fermi- 
Dirac probability distribution using zu m =7.5 identical to 
the experimental data limits for purposes of comparison, and 
the third row contains the partial moments evaluated for the 
probability distribution extracted from the experimental data 
also using zu m = 7.5. Very good agreement is seen. 

On adopting this procedure of multiple moment com- 
parison for — z\[ m < z < zii m , we find optimised values of 




FIG. 4: a Data (circles) from Fig. [3] plotted versus 77 = 
27rkBTm* /eB, where m* =1.676 ra e . The data from the lower 
inset to Fig. [3] is included in this plot after renormalizing 
its amplitude relative to the contactless conductivity data at 
77 « 0.15. The red line corresponds to a linear interpolation, 
b Numerical Fourier transform of a (77) in (a) (red line) ob- 
tained using Eqn.[T] plotted over a restricted range of positive 
z — e/kBT, with the Fermi-Dirac (black), Gaussian (green) 
and Lorentzian (blue) distributions shown for comparison as 
dotted lines. 



ao and m* = 1.676 ± 0.001 for which the second (vari- 
ance) and fourth (kurtosis) moments of the experimen- 
tally obtained probability distribution concides to 4 and 
3 significant digits with the Fermi-Dirac probability dis- 
tribution, respectively. Good agreement continues to be 
found with the higher order moments, which we list in 
Table I up to /iio- In Fig. [5] we compare the Fermi-Dirac 
distribution and the illustrative Gaussian and Lorentzian 
probability distributions against that obtained experi- 
mentally. This exercise demonstrates the importance of a 
wide temperature range in testing the goodness of fit of a 
probability distribution model— we see that the deviation 
from the Gaussian and Lorentzian models only appears 
above z = 2. The deviation from an ideal Fermi-Dirac 
distribution occurs only above z — 6, but falls within the 
experimental uncertainty. 

Our study therefore indicates a high degree of corre- 
spondence between the elementary excitations in under- 
doped YBa2Cu306+^ (x ~ 0.56) and those of fermionic 
quasiparticles governed by a Fermi-Dirac distribution. 
Hence, long-lived, robust fermionic quasiparticles ex- 
ist over a broad range of temperatures (100 mK < 
T < 18 K). At first glance one may find this result 
surprising, especially in a system such as underdoped 
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FIG. 5: a Comparison of the fit probability distributions 
shown in Fig. [4] with that obtained from experiment. The 
Fermi-Dirac (red line), Gaussian (green dotted line) and 
Lorentzian (blue dotted line) probability distributions are 
shown on the y-axis as a function of the experimentally ob- 
stained probability distribution |f'(;z)l on the x-axis. The cir- 
cles locate integer values of z (the black dotted line is a guide 
to the eye) , with the error bars indicating the extent to which 
experimental uncertainties (chiefly associated with the sample 
temperature) can cause |f'(z)| to depart from |fpD(X)l- 

YBa2Cu306+x which is close to the Mott insulating 
state. We relook at whether we indeed have reason to 
find this result unexpected, and investigate possible rea- 
sons why Landau Fermi liquid behaviour prevails. We 
consider alternative scenarios that may be applicable to 
the ground state of underdoped YBa2Cu306+^, and in- 
vestigate whether they would be consistent with our ex- 
perimental observation of Fermi-Dirac statistics. In ad- 
dition to non Fermi-Dirac statistics, the temperature de- 
pendence of the oscillation amplitude could deviate from 
that shown in Fig. [T] due to a non-constant value of go. 
An energy gap in the electronic structure near the Fermi 
energy of magnitude comparable to hco c could yield a 
non-constant value of go, such as in the case of symmetry- 
breaking groundstates (e.g. density wave ordered states 
with Fermi surface reconstruction due to translational 
symmetry breaking and the superconducting vortex liq- 
uid) . Any additional temperature dependence of the am- 
plitude arising from these effects, however, appears not to 
be large enough to lead to a detectable deviation in our 
experiment from the temperature-dependent amplitude 
expected from Fermi-Dirac statistics. 

Various other models would not be expected to lead 
to Fermi-Dirac statistics - a notable example being 
fractional quantum Hall effect in very high magnetic 
fields. In the case of underdoped YE^CuaOe+x, mod- 
els yielding oscillations unrelated to Landau quantisa- 



tion (e.g. Ref. [lf|) would not be expected to result in a 
temperature dependent oscillation amplitude that agrees 
with Fermi-Dirac statistics, and are therefore inconsis- 
tent with our experimental findings. To consider the 
applicability of various other non Fermi liquid models 
proposed for YBa2Cu306+^ (e.g. Ref. [8]), theoretical 
predictions are required for the energy scale at which a 
deviation from Fermi-Dirac statistics is expected. Such 
ground states could potentially go undetected in the 
current study if the deviation from Fermi-Dirac statis- 
tics takes place over an energy scale lower than that 
~ 100 mK of the lowest accessed temperature. Similar 
considerations may apply to the superconducting vortex 
liquid groundstate in which regime the current exper- 
iments are performed. It will be of interest to deter- 
mine conditions under which a breakdown of Fermi-Dirac 
statistics could occur for this groundstate. 
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